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ABSTRACT 


This  work  covers  two  major  topics.  Firstly,  crystal  growth  experiments  on 
a-HgS  are  reported.  The  aim  here  was  to  use  the  THM  method  of  solution  growth  to 
grow  large  crystals  from  solvents.  So  far,  no  large  crystals  have  been  obtained 
but  the  general  utility  of  as  a  solvent  for  the  growth  of  smaller  cinnabar  crystals 
could  be  shown.  Also,  proustite,  lead  cesium  chloride  and  bismuth  sulfide  were 
grown  using  the  THM  method. 

Secondly,  a  method  for  the  calculation  of  electrooptic  coefficients  of  materials 
using  only  basic  parameters  was  developed.  Phillips  electronegativities,  crystallo¬ 
graphic  data  and  electronic  polarizabilities  of  the  atoms  are  used  as  inputs  in  a  com¬ 
puter  program  which  is  described.  The  program  outputs  the  elcctrooptic  coefficients 
of  cubic  and  noncubic  structures.  The  major  test  cases  were  compounds  of  the  zinc- 
blende  structure  and  extremely  good  agreement  between  calculated  and  experimen¬ 
tally  determined  eleciroopdc  coefficients  is  obtained. 
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I.  INTRODUCTION 

Recently  it  could  be  shown,  in  a  number  of  instances,  that  certain  methods 
of  solution  growth  are  capable  of  growing  materials  of  interest  for  electrooptics 
which  are  otherwise  difficult  to  prepare.  The  methods  of  solution  growth  used  are 
the  "Traveling  Solvent  Method"  (TSM)  and  the  "Traveling  Heater  Method"  (THM), 

the  latter  one  particularly  applicable  to  the  growth  of  larger  crystals.  Thus,  growth 

1  2  3 

of  ZnO,  (PbxSrj_x)TiOj,  and  CdTe  was  accomplished  and  for  the  latter  two  cases 
electrooptic  modulation  properties  were  demonstrated. 

Since  solution  growth  methods  are  often  very  powerful  for  the  preparation  of 
complicated  compounds,  in  particular  if  solid  state  transformations,  high  vapor 
pressure,  and  deviations  from  stoichiometry  are  involved,  it  seems  reasonable  to 
test  THM  in  particular  as  a  method  for  growth  of  various  complex  sulfides  and  other 
compounds  which  are  of  interest  as  potentially  useful  modulator  materials.  .In  par¬ 
ticular,  the  growth  of  sizable  crystals  of  cinnabar  (arHgS)  has  always  presented 
some  problems  because  of  its  phase  transformation  at  344  °C  into  metacinnabar,  the 
cubic  form.  Thus,  growth  of  a-HgS  crystals  with  sizes  sufficient  for  electrooptic 
usage  was  the  main  objective  of  the  program. 

Although  in  some  cases  a  unique  necessity  for  crystal  growth  from  solution 
exists,  this  is  generally  not  the  case.  Alternative  methods  of  single  crystal  prep¬ 
aration  often  exist.  Nevertheless  growth  from  solution,  if  properly  conducted  such 
as  in  the  THM  arrangement,  leads  often  to  far  higher  perfection  and  purity  in  the 
final  crystal.  This  is  due  to  the  fact  that:  (1)  growth  occurs  slowly  under  near 
equilibrium  conditions,  and  (2)  impurity  segregation  is  often  facilitated  if  a  solvent 
is  used. 

Apart  from  these  experimental  objectives,  a  major  theoretical  objective  has 
been  adopted.  This  concerns  the  absence  of  reliable  models  for  the  calculation  of 
Various  electrooptic  properties  of  materials.  Thus,  work  on  a  simple  model, 
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based  on  readily  accessible  Information  of  a  physiochemical  and  crystallographic 
nature,  which  allows  reliable  calculation  of  electrooptic  coefficients  and  other 
optical  properties  of  crystalline  materials  has  been  carried  out. 


II.  WORK  IN  THE  FIRST  SIX  MONTHS  OF  THE  CONTRACT 

This  work  is  described  in  detail  in  Seminanual  Report  No.  1  on  the  present 
contract,  which  was  issued  in  January  1971  and  will  only  be  summarized  here. 

The  major  amount  of  experimental  work  in  the  first  half  of  the  contract  was 
directed  towards  finding  a  suitable  solvent  for  the  THM  growth  of  cinnabar.  In  the 
course  of  these  experiments,  a  compound  HggSglg  was  found  and  the  phase  diagram 
of  the  HgS-Hglj  system  was  delineated. 

Nevertheless,  the  only  suitable  solvents  found  for  a-HgS  were  polysulfides 
of  sodium  and  potassium  and  Initial  growth  experiments  using  these  solvents  were 
carried  out  without  too  much  success. 

Also  PbCsClg  was  grown  from  PbClg  melts.  Although  no  single  crystals 
resulted,  the  regrowth  of  polycrystalline  boules  clearly  demonstrated  that  further 
cultivation  of  the  method  using  a  PbClg  solvent  would  undoubtedly  lead  to  the  growth 
of  large  crystals.  Since,  however,  the  material  does  not  seem  to  be  of  too  much 
interest  at  present,  further  experimentation  was  discontinued. 

Also,  several  experiments  on  the  growth  of  proustite  from  ASgSg  rich  solu¬ 
tions  were  undertaken.  These  experiments  were  quite  successful  in  producing  siz¬ 
able  clear  crystals  in  several  instances.  There  were  indications  from  optical  trans¬ 
mission  measurements  that  through  the  use  of  an  ASgSg  rich  solvent  the  absorption 
peak  around  10  m  might  be  minimized  or  avoided.  The  disappearance  of  that  absorp¬ 
tion  was  not  reproducibly  observed  in  all  samples  grown,  but  one  might  be  led  to  the 
conclusion  that  it  could  be  caused  by  a  native  defect. 

Further  work  in  the  first  8-month  period  of  die  contract  was  concerned  with 
a  literature  search  which  covered  a-HgS,  proustite,  lead-cesium  chloride,  andother 
ternary  noncentrosym  metric  compounds  which  might  be  of  interest  in  electrooptics. 

Finally,  a  major  amount  of  work  was  devoted  to  the  development  of  a  theory 
which  allows  calculation  of  die  electrooptic  coefficients  from  basic  parameters.  The 
method  usesas  basic  information  only  crystal  structure  data,  lattice  polarizabilities 
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and  electronegativity  differences.  Computing  the  electrooptic  coefficients  for  the 
cubic  zincblende  structure,  very  good  agreement  between  the  computed  values  and 
measured  values  was  indeed  obtained.  Extension  of  the  method  to  noncubic  struc¬ 
tures  was  also  started  and  is  described  in  detail  in  the  body  of  this  report. 
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and  tj j  and  arc  appropriately  defined  lattice  sums.  This  expression  has  been 
derived  for  a  rigid  but  polarizable  ionic  model  under  the  assumption  that  the  internal 
and  external  fields  are  parallel  which  is  true  for  cubic  materials  or  special  cases  in 
other  symmetries.  For  example,  in  a  hexagonal  crystal,  when  the  applied  electric 
field  is  parallel  to  the  e-axis  or  in  the  basal  plane,  the  internal  field  will  be  in  the 
same  direction.  An  isotropic  dielectric  constant  has  been  tacitly  assumed. 

In  the  more  general  ease,  an  expression  for  die  internal  fields  at  the  lattice 
sites  has  been  derived.  Without  going  into  details,  it  can  be  shown  that: 


j  (t+ai)  E  + 


E 

j,  n,  m,  ( 


.  (X,+xpx:  €„V  r 

*J  flJ  ' 


5<V<Yj 


r*  F* 

rlj  h) 


t..  ■  2  (f„  ■  ?,)  J  * 


1J 


ij  1 )'  rIJ  FJ 


<«> 


No  simple  general  expression  for  the  dielectric  constant  such  ac  Eij.  (1)  can  be 
written  down.  Note  that  here  the  dielectric  cons  taw  is  treated  as  a  tensor  and  Xj 
probably  should,  bu  tas  not. 

Returning  to  tl*?  more  restricted  case,  the  first  term  on  the  right  hand  side 

in  Eq.  (2)  for  r4,  is  th<  contribution  to  the  dielectric  constant  by  the  polarizability  of 
*J  2 

the  ions.  For  optical  frequencies  where  c  ■  n  ,  the  ionic  part  of  Xj  will  dlsppcnr 
(i.e.,  Xj  -  0) .  The  second  term  is  the  Lorentz  correction  which  can  contribute  to  the 
dielectric  constant  if  the  structure  of  the  lattice  allows  a  large  value  for  t^  which 
occurs  for  ferroelectric  materials  such  as  BaTiOj.  The  third  term  produces  an 
electric  field  dependence  of  the  dielectric  constant  or  at  optical  frequencies  gives 
rise  to  a  linear  electrooptic  effect.  This  term  vanishes  for  centrosymmetric  ma¬ 
terials.  The  next  higher  order  term  (which  has  not  been  written  out)  can  produce 
the  dielectric  constant  nonlinearity*  and  quadratic  electrooptic  effect  observea  in 
materials  such  as  SrTiOj.® 

Let  us  now  consider  the  specific  case  of  the  wurtzite  structure.  Table  i  shows 
the  location  of  the  four  atoms  in  the  hexagonal  unit  cell.  Also  shown  are  the  vectors, 


Table  I.  Coordinates  of  Atoms  and  Lattice  Vectors 
In  the  Wurtzite  Hexagrnal  Unit  Cell  -  AB 


Coordinates  of  Atoms 


Atomic  Position 

X 

y 

z 

A-atonu  1 

ma 

na 

f  c 

9 

(m  +  J)  a 

(n  +|)  a 

(1  4  j)  C 

B-atoms  9 

ma 

na 

(f  +u)  c 

4 

(m  +  J)  a 

(n  +  j)  a 

(1  +  J  +u) 

Coordinates  of  Lattice  Vector 


Lattice  Vectors 

X 

y 

z 

rll  *  r9S  *  r99  ■  r44 

ma 

na 

Ic 

rII  -  r94 

(m  +J)  a 

(n  *  f)  * 

U  c 

r19  *  r84 

ma 

na 

(f  +u)  c 

rl4 

(m  +  j)  a 

<"  -  j>  • 

(I  +u  -  j)  c 

rI9 

(m*j)a 

(n  +  j)  ■ 

(I  4U  -  J)  C 

In  an  Ideal  wurtzite  lattice,  c/a 
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r^,  that  connect  the  different  lattice  sites.  Table  II  shows  the  lattice  sites  for  cinna¬ 
bar  using  the  three  molecule  unit  cell.  (Nothing  further  has  been  done  with  calcula¬ 
tions  on  cinnabar  this  report  period,  and  we  will  restrict  ourselves  to  wurtzitc.) 
Writing  out  Eq.  (1) : 


c  -  1 


I 

5 


2  <XA  +  Xg)  +  (2tjj  4  2tw  +  2tjj  4  tJ4 


+ 128^  +  ^  +  ®8|1  +  *®I2  **s13  +  S14 


+  822)^a<XA  +  XA)  +^B(XB+Xil 


where 


ij 


TT  S 

1  n,  m,l 


ru 


•u 


ir  E 

B  n,  m,  f 


3(fi-fn)  (8<n.f(|)J  -Sr*] 

rv 

U 


(5) 

(«) 

(7) 


where H  is  either  along  the  z  (optic)  axis  or  in  the  basal  plane.  Use  has  been  made 
of  die  symmetry  of  t^  and  s^  to  reduce  the  number  of  independent  components.  The 
solution  for  an  arbitrary  direction  is  not  obvious  as  Eq.  (4)  would  have  to  be  used. 
Also  die  anisotropy  of  c  has  been  neglected. 

We  will  choose  as  the  coordinate  system  that  of  die  axes  of  the  hexagonal  unit 
cell  which  is  nonorthogonal.  Thus: 

*2  2  2  2 

r  •  x4  +  y*  •  xy  +  z 

if  x,  y,  and  z  are  die  diree  components  off. 

Also: 

r*fl  *  xnj  4  y -  y  (xn^ynj)  4  ziij 
where  n^,  n^,  and  n^  are  the  three  components  of  n. 


-  • 


9 


Table  11.  Coordinates  of  Atoms  in  Cinnabar  Using  A  Trimolecider 
Hexagonal  Unit  Cell 


Coordinates  of  Atoms 


Atomic  Position 

X 

y 

2 

Hg-atoms  1 

(m  +  u)  a 

n 

«  +j)  * 

S 

ma 

(n  4  u)  a 

(1  +  §)  c 

3 

(m  -  u)  a 

(n  -u)  a 

fc 

S-atoms  4 

(m  4  v)  a 

na 

(1  +J)  c 

5 

ma 

(n  4  v)  a 

(1  +$)  c 

6 

(m  -  v)  a 

(n  -  v)  a 

(1  +j)  C 

u  -  0.710  a  -  4.140 

v  -  0.406  c  -  0.406 
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The  computer  program  to  perform  these  sums  Is  shown  below. 

Generally  the  calculations  of  lattice  sums  can  be  performed  in  a  straightforward 
manner.  Using  a  computer,  the  calculations  are  made  over  successively  larger  and 
larger  volumes  about  the  origin  which  allows  using  a  Neville  plot  to  extrapolation  to 
infinite  volumes.  A  Neville  plot  is  plot  of  the  sum  versus  N~*  where  N  is  the  total 
number  of  atoms  over  which  the  sum  is  performed,  thus  extrapolating  to  where  N~* 

7 

is  zero  gives  the  limit  of  the  sum.  This  approach  is  particularly  effective  when  the 
terms  have  the  form  of  r"n  where  n  >  S.  If  n  s  3,  convergence  can  be  a  problem  un¬ 
less  some  care  is  taken  to  arrange  the  terms  so  that  contributions  from  negative  and 
positive  terms  nearly  cancel.  In  this  program  the  use  of  a  spherical  volume 
worked  fairly  well.  This  was  checked  since  it  is  known  that  (2tjj  +  2tJ2  +  2tjj 
4  tjj  4  tjj)  should  be  zero  for  wurtzite,  and  when  the  calculation  was  made  this 
quantity  was  small. 

The  first  part  of  the  program  inputs  the  value  of  c/a,  u,  and  the  largest  value 
of  r  over  which  the  sum  is  to  be  made.  The  four  separate  sums  are  calculated  at  the 
same  time  and  are  identified  such  that  the  index  1  corresponds  to  11,  2  to  12,  3  to  13, 
and  4  to  14.  The  sum  is  then  made  over  each  atom  by  calling  die  subroutine  appro¬ 
priate  to  the  external  electric  field  being  applied  parallel  or  perpendicular  to  die 
c-axis.  Typical  output  is  shown  where  the  sum  is  shown  for  each  term  over 
larger  unit  cells.  In  this  particular  case  for  an  electric  field  along  the  c-axis, 
Sj2  and  Sjj  are  zero.  Eq.  (1)  or  (8)  can  be  used  to  calculate  die  change  in  dielectric 
constant  or  electrooptic  effect  since  n*  ■  c  and  Xj  -  0. 

2.  Computer  program  for  lattice  sums 

The  following  program  was  used  to  calculate  the  lattice  sums. 

r  psnrpvi  to  rai/M.^TF  l.aTTlCF  For  a  "Mr*  J tf  **t- 

n j  "rm-1  'u*) 

Tvp-  1  -» 

I*  ///'  TVPVT  r/a, 

rr 

r.  r.n*  I"  t*’~  ’’itjo  nr  n/.a  <t  i.rva. 

C  I'Y  Tf  T,,r  n,"^Tur’'*T  TP*’  0VFn  v»'7C”  T'T  **'l  J  !<■  T n  ***  Mr,# 

C  M  tV ’  •'T^rpir-  pa^'MFTF**  fll’P^PyT  MTH.v  FP'fM.  Tp 
F*F:na*r*nr  m*** 
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a  n  Avr  r  Ar“r*  T»*r  «'ATPr  p f  a  to****  i*  T""  r*>Jw-TTT'.r  rr? 1 

*(*>rl./3. 

A(3>r". 

A(A): A<?) 

PCIJrP. 
n(?)s- V?) 

n(A) Tr f?) 

r.(  i  )  =  «. 

C^Jr.? 

C(3)r!i 
C(A)rU-.f 
Kr * 

TyPF  jr  F'n,M,rn!\f-‘Ay#K 
3*  FnR1AT(/3Fir.P,'»Tf  //) 

Tvpf  3? 

33  rOR*|AT  ( //  *  TFPi  TT  c*r  A  n 

1  r  TOM**  *  ) 

7"  KrK-H 
XKrVF+1 
TTTrF. 

rn  j*  !ri,A 

TT=?. 

rr-*\ 

NinwrrF 

AArAU) 

',0SP(J, 

criscm 

rn  af 

yfsMN-l-F 

no  af  'nsi  ** 

X’l:  M'l- 1  -v 

no  af  u.sI.fk 

X1.rU.-l-F 

r  thf  coa’trjphtj  n\'  j*  FAi.cui.4TFr  fo*>  fac»»  tfr  »  r»  a 

C  r PM FRF  OF  OAPT'r  V , 

cai.i.  TF*n<x*’,y*i,yi..,M,p«,nn,?,'*> 

TTrTT*T 
FFr  Ff+r 
AF  CftMT!  fj||F 
TTrTT*FAC 
fCrfeaFAC/P. 

TTTrTTT4-TT 

V  TVPC  r*,,T,TT,,,r, AA,pp,rnf •.'inf**' 

F"  F0R1AT(3y. I3t?F I  A.* ,3FP,A  f I*  ) 

TVPF  PP.TTT 

rf  FnpM4T(/,?y,  jam  <*mi  if  tt  tr  ) 

TVPF  F? 

F«  FPM'IAT  (/) 

IF(F.I.T.'TAy)  R0  TP  7F 
rn^TIfic 

rvn 
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The  subroutine  for  the  case  the  electric  field  is  along  the  c-axis  is  the 
following: 


r  cai.cmi.ati  nr  Tr.p,M*  in  the  cm  por  c-Ayic  '*MiR7iTr. 

*,iins»nirT!fi"  Trm<y,v,7,  A,p,c,T,r> 
cptwv  c.A,yt"»,v 

c-  n 
•  *  •  • 

T  z  ?. 

I  F(x*y-»-v*v'-*.7'‘r7*r:A*CA-y*v,  ni  ,vk?  >pET'ipfj 

pr>p.  (yj./\)*(y+ft)4.(Y-*.R)*(V+P)-(X-*-A)*(V+n)4-(7+r>*(7+r)*rA*P.A 

I rcRPn.rn.P. )  °rTMR*.t 
Vs 

■? n v -  (7+P  )+C  A 

t=  (RrE-5.*pr*'*n-‘’)/ri'r**-»* 

r=  3 . +p  n  M*  ( ■> .  *°  r  ,,*p  n  n -3 .  *p  rp  >  /p  rp  ♦*  3  • * 

RFTI'RV 

r\r 


Typical  output  is  shown  below: 


ivpilT  M,  C/A,  *1A v f  .373  I.C3PCC3I  f 


P.1IP5AP  P.373PPP  |,^3?Rp3  K  p 


TEP'l  TT 

.  A  «* 

C 

I  nv«* 

I  r\*7^At>rv 

P.PPPPPF+PP 

n ,  ppm  p  pppp 

P. PPPP 

r 

9.  -P.3?.3irp+np 

P.39435E-PI 

ft. 3333  -  P  !  3  3  3 3 

ni  ^  ppp 

7 

3  -p.opAPcr+p'' 

P. 1CA3PF+P1 

P.PPPP  P.PPPP 

P  *  3  7?  ^ 

7 

A  P.I!*33S>PI 

p.3PArpr+?i| 

"•.3333  -P.3333 

- fl.  | P* P 

7 

CH  07  TT  Ic 

P.C1337F+PP 

1 

p.??c«p*>pp 

P.PPPPPF+PP 

P.PPPP 

P.PPPP 

3° 

v 

-P.K/3CPE+P" 

-P.I'MACF+PP 

'’.3333 

-7.3333 

*» 

33 

3 

-P.  !  IcmC+PI 

7.1  AP5E+PI 

o  mnj  j»n 

P, PPPP. 

P..375P 

33 

A 

P. ! I ! PFE+Pl 

P.3  ICC f r+P| 

'•13333 

-P.3333 

-P,  imp 

3.3 

cm  OF  TT  T r 

-P. 1757AE+PP 
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1 

0,3  90Q  |  T>00 

O.OOOOOF+OO 

0.0000 

ft  ft  ft  OJ  ft 

rt  onirtp 

74 

2 

-0.55F57F+PP 

-0.4227FF-02 

0.3333 

-oi3SS3 

0)5000 

75 

3 

-0.  1  lOOFF-t-OI 

0. 1 7002F+01 

0,0000 

0, 0000 

0.3750 

75 

n 

O.l  I402.F+0! 

0.317FIF+0I 

0.3333 

-0.3333 

-0. |?50 

73 

F'll  OF  TT  IF 

0. *07**2-01 

1 

0.577*5F+00 

-0. 13I01F-0? 

CT<  '»!*'»(? 

*  (*  ft 

Of»  17 

2 

-0.574O3F+00 

-0.4OJ  74F-0?. 

0.3333 

-0)3333 

0.5000 

209 

3 

-0. 1 M53F+OI 

0.  ISO 7?.F-*-01 

0, 0000 

0.0000 

0.3750 

O 

4 

O.l  !!7i»F^I 

0,31731 F+0| 

0.3333 

-0.3333 

-0. 1250 

209 

FIM  OF  TT  !*■ 

-0.25IO7F-0I 

1 

0.5  7FFIF4.00 

-0.25*752-0* 

0, 0000 

0, 0000 

0. 0000 

3«?4 

<y 

-0.5707FF+00 

0.FF4A4F-03 

0.3333 

-0.3333 

0.5000 

3*3 

3 

-0.1 1 A ! 7F+0 1 

0. I7017F+01 

0,0000 

(T»  finO"1 

0, 3750 

303 

4 

0.1  1  I7FF+0I 

0.3l7|FF-t-0| 

0.3333 

-0)3333 

-0. 1250 

30  3 

FIM  OF  TT  I* 

-0,1*3032-01 

1 

O.5fi«>,»7F4-0P 

0.  |*«|*F-0* 

0, 0000 

pi  nnnni 

0.0000 

*3* 

2 

-0.5*0 !? F+nn 

-0.2051 1F-03 
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-0)3333 

0.5000 

r  37 

3 

-o.l  150*04.0! 

0. 1  70nof>0! 

0.0000 

n 

0.3750 

*37 

4 

0.1  I  0*  1  r>0  1 

0.3I722F+0I 

.0,3333 

-0)3333 

-0. 1250 

*37 

PIH  OF  TT  J«?  -O^IFFOF-O! 
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3.  Crystal  growth  by  THM:  A  detailed  qualitative 
process  analysis* 

Any  analysis  of  THM  must  consider  the  fact  that  the  diffusion  process 
is  convection  aided,  as  well  as  external  experimental  variables  such  as  heat  flow 
through  solvent  and  the  adjacent  solids,  etc. 

Recently,  such  an  initial  qualitative  analysis  was  carried  out  by  Tyco  per¬ 
sonnel  working  under  several  contracts  which  use  the  THM  method  (including  the  one 
on  which  we  report)  with  the  help  of  Professor  Bruce  Chalmers  of  Harvard  University, 
who  is  a  consultant  to  Tyco  on  crystal  growth  matters.  The  results  of  this  initial 
analysis  are  reported  in  this  section. 

Essentially  what  is  being  done  is  to  analyze  the  conditions  for  steady  state 
THM  growth  under  various  experimental  constraints.  Without  serious  loss  of  pre¬ 
cision,  the  following  assumptions  are  made: 

1.  The  relevant  liquldus  line  is  straight. 

2.  The  solidus  line  is  vertical. 

3.  Each  solid-liquid  interface  is  at  the  liquidus  temperature  for  the 
liquid  at  the  interface.  (This  implies  that  growth  is  slow  enough  for  kinetic  super¬ 
cooling  to  be  neglected.) 

The  following  sequence  Is  adopted: 

1.  Description  of  the  process 

2.  Solute  and  temperature  distribution  in  the  liquid  zone:  (a)  when 
transport  is  by  diffusion  alone,  (b)  when  stirring  is  completely 
effective,  and  (c)  when  stirring  is  partially  effective 

3.  Effect  of  temperature  fluctuations 

4.  The  shape  of  the  Interface. 

The  Process 

For  the  present  purpose,  the  process  can  be  represented  by  Fig.  1 
in  which  a  liquid  zone  moves  upwards  relative  to  the  growing  crystal  C.  The  relevant 
part  of  the  phase  diagram  is  shown  in  the  right  hand  part  of  the  figure.  The  liquid 

*F.  Wald. 
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/one  AB  has  a  thickness  Z,  which  is  essentially  independent  of  temperature  and  speed. 
The  liquid  at  A  has  a  composition  Cj,  and  that  at  B,  Cj.  There  is,  therefore,  a  con¬ 
centration  gradient  of  average  value  (Cj  -  Cj)/Z  through  which  transport  takes  place, 
tending  to  increase  the  concentration  at  B  and  to  reduce  it  at  A.  This  tendency  is  pre¬ 
vented  by  crystal  growth  at  B  and  solution  at  A.  The  process  will  continue  so  long  as 
Tj  >  Tj.  This  temperature  difference  can  be  maintained  in  two  ways:  (1)  by  im¬ 
posing  a  temperature  gradient  on  the  system,  and  moving  the  temperature  profile 
relative  to  the  sample  at  an  appropriate  rate,  or  (2)  by  supplying  heat  assymmetri- 
callv  to  the  liquid  zone  by  continuous  relative  motion  and  extracting  it  from  the  two 
solid  regions.  The  heat  source  is  moved  relative  to  the  sample.  (This  latter  is  truly 
the  essence  of  the  THM  process.) 

It  can  be  shown  that  the  second  method  (heated  zone)  offers  much  more  flexi¬ 
bility  than  die  first  (moving  gradient) .  Assuming  (for  the  moment)  that  transport  is 
purely  by  diffusion,  the  following  conditions  would  exist.  With  a  fixed  temperature 
gradient  [Case  (1)  ]  and  a  fixed  zone  length,  Tj  -  Tj,  Cj  -  Cj  and  Z  would  be  fixed; 
this  would  produce  a  flux  [D  (Cj  -  Cj)  /Z]  through  the  zone,  which  would  define  the 
rate  of  growth.  Unless  the  temperature  gradient  is  moved  at  the  same  speed,  the 
solvent  zone  would  move  continuously  up  or  down  the  temperature  gradient,  and  no 
steady  state  would  be  possible.  In  the  second  case,  on  die  other  hand  (Fig.  2.),  a 
fixed  zone  length  would  be  compatible  with  a  wide  range  of  temperature  differences 
(from  zero  to  T j  -  in  Fig.  2),  and  a  correspondingly  wide  range  of  values  of 
Cj  -  Cj.  The  flux  through  die  zone  can,  therefore,  adjust  itself,  by  suitable  position¬ 
ing  of  die  zone,  to  cause  crystal  growth  at  any  speed  within  a  wide  range.  There  is, 
therefore,  considerable  latitude  in  the  speed  at  which  continuous  steady  state  growth 
can  occur  for  a  given  temperature  profile.  This  will  change  only  in  detail  but  not  in 
significance  if  the  difference  in  thermal  conductivities  between  the  solvent  and  die 
solid  is  to  be  taken  into  account. 

Solute  and  Temperature  Distribution  in  the  Liquid  Zone 

Three  types  of  defects  may  occur:  (1)  "stray"  crystals,  (2)  in¬ 
clusions,  and  (3)  cellular  substructures.  These  can  all  arise  from  cither  general  or 
local  supercooling  of  the  liquid  near  the  growth  interface.  General  supercooling  is 
that  which  occurs  right  across  the  section,  as  opposed  to  local,  which  does  not.  Since 
heat  is  applied  to  the  liquid,  the  latter  must  (if  we  assume  a  planar  interface)  be  at 
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Pig.  S.  Consequences  of  using  an  asymmetrical  temperature  distribution 
in  the  solvent  as  in  the  THM  method 
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a  higher  temperature  than  TjJ  the  only  supercooling  that  can  occur,  therefore,  would 
be  constitutional.  We  will  consider  the  relationship  between  composition  and  tempera¬ 
ture  for  a  planar  interface,  making  various  assumptions  about  the  mechanism  of 
transport.  Three  assumptions  will  be  made:  (a)  transport  by  diffusion  alone  (no 
mixing, )  (b)  completely  effective  mixing,  (c)  complete  mixing  except  for  a  stagnant 
diffusion  zone  at  each  end. 

a.  Diffusion  alone:  Steady  state  diffusion  from  interface  1  to 
interface  II  requires  a  uniform  gradient  of  composition.  If  the  temperature  difference 
is  maintained  by  a  uniform  temperature  gradient  (Fig.  9),  the  temperature  of  the 
liquid  would  be  equal  to  its  liquidu3  at  all  points.  If  the  heat  is  supplied  to  the  liquid, 
the  distribution  of  temperatures  [actual  (Ta)  and  liquidus  (T^)]  is  shown  in  Fig.  4. 

The  actual  temperature  is  everywhere  above  the  liquidus,  and  there  is  no  supercooling. 

b.  Complete  mixing:  This  would  give  a  composition  profile  as 
shown  in  Fig.  5.  The  liquidus  temperature  would  be  as  shown  in  the  center  portion 
of  the  figure  which  also  shows  a  uniform  gradient  (unrealistic  in  a  folly  mixed  liquid) 
and  a  temperature  profile  TA  that  might  exist  in  a  stirred  heated  liquid.  In  the  latter 
case,  there  is  a  small  supercooled  zone. 

c.  Partial  mixing:  The  temperature  distribution  for  the  heated 
liquid  case  would  be  as  shown  (Fig.  6) ,  the  stagnant  zones  having  high  temperature 
gradients.  There  could  be  a  small  supercooled  zone,  but  the  fact  that  the  "mixed 
region"  applied  to  temperature  as  well  as  to  solute  suggests  that  a  substantial  region 
of  supercooling  is  not  to  be  anticipated. 

The  general  conclusion  of  this  section  is  that  under  realistic  experimental 
conditions  and  assuming  a  planar  interface,  there  is  no  severe  constitutional  super¬ 
cooling  in  THM. 

Effect  of  Temperature  Fluctuations 

We  will  consider  the  possible  effect  of  temperature  fluctuations 

that  take  place  too  fast  for  diffusion  to  occur  over  distances  comparable  with  the 

-B  1  2 

thickness  of  the  liquid  zone  (if  Z  *  1  cm,  D  ■  10  cm  /sec;  then,  since  d  ■  OT, 
t  *  10®  sec) .  Therefore,  fluctuation  of  period,  say  1  minute  or  1  hour  would  qualify; 
but  a  period  of  1  day  would  not.  Assume  that  steady  state  has  been  achieved  ("diffusion 
only"  case);  the  liquid  at  1  is  at  Cj,  etc.  (Fig.  7).  Now  lower  the  temperature  to  Tj' 
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CONCENTRATION  OF  UQUtO  UOUtOUS  TEMPERATURE  — » 

Pig.  S.  Temperature  and  concentration  In  the  solvent  zone  for  the 
equilibrium  case  (pure  diffusion) 

xl 

T. 

Ta 

I  ^ 

TEKII?— » 

Pig.  4.  Relation  of  liquidus  temperature  and  actual  temperature  for  the 
case  of  pure  diffusion  in  THM 
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Pig.  5.  Relation  of  concentration  and  temperature  for  the  case  of  a  fully 
stirred  liquid 


Pig.  6.  Temperature  and  concentration  distribution  for  the  case  of  a  partially 
stirred  liquid 
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and  Tj*.  The  liquid  close  to  die  interfaces  will  begin  to  change  towards  Cj’  and  Cj'; 
after  a  short  time,  die  composition  profile  would  be  curve  B  in  Fig.  7;  the  liquidus 
temperatures  will  be  shown  on  the  left  in  Fig.  8.  The  actual  temperature,  however, 
will  be  as  shown  on  the  right;  a  supercooled  region  can  exist  near  interface  11  if  the 
temperature  change  is  relatively  fast.  This  effect  is  reduced  by  a  steeper  tempera¬ 
ture  gradient  in  the  liquid,  but  this  gradient  is  limited  by  the  (usually  high)  thermal 
conductivity  of  the  solvent. 

Interface  Shape 

The  foregoing  discussion  suggests  that,  if  the  underlying  assump¬ 
tion  of  planar  interfaces  is  valid,  substantial  supercooling  is  unlikely  but  not  im¬ 
possible.  However,  the  interface  may  not  be  planar.  We  will  consider  first  the  shapes 
of  the  isothermal  surfaces  that  characterize  the  cases  (a)  when  heat  transfer  from 
source  to  sink  is  by  conduction  alone,  and  (b)  when  there  is  sufficient  mixing  that  the 
liquid  is  isothermal  except  in  die  two  boundary  layers,  where  heat  transfer  is  by 
conduction.  This  will  be  followed  by  a  discussion  of  the  conditions  under  which  the 
interface  may  deviate  from  the  isotherm,  and  the  possible  consequences. 

a.  Heat  transfer  by  conduction  alone;  In  the  (unrealistic)  case  in 
which  the  solid  and  die  liquid  have  equal  thermal  conductivities,  the  thermal  flow  lines 
and  the  isotherms  are  shown  in  Fig.  9-a  if  the  heat  source  is  shorter  (along  the 
sample)  than  the  heat  sink,  9-b  if  the  source  is  longer,  and  9-c  if  they  are  equal. 
Cases  b  and  c  can  be  achieved  in  detail  only  by  adding  other  flow  lines;  but  the  im¬ 
portant  feature  of  b,  the  existence  of  isotherms  that  are  convex  upwards,  seems  to 
be  possible  in  a  system  like  the  one  under  consideration.  It  will  be  noted  that  Fig.  9-c 
includes  the  only  case  of  a  planar  Isotherm.  It  is  not  clear  whether  this  is  possible 
in  die  actual  system,  because  die  required  symmetry  requires  isotherms  such  as  AB, 
which  implies  another  heat  source  lower  down. 

The  assumption  of  equal  thermal  conductivities  is  unrealistic.  In  general  the 
ratio  is  much  larger  than  1,  which  seems  to  mean  that  conditions  close  to  Fig.  9-c 
would  require  the  heat  source  to  be  four  times  as  long  as  the  heat  sink.  This  suggests 
that  conditions  like  9-b  can  occur,  and  therefore  that  convex  isotherms  are  possible; 
they  are  not  under  the  conditions  of  Fig.  9-a,  which,  without  the  conductivity  ratio, 
would  be  like  the  actual  case. 
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Pig.  7.  Effect  of  temperature  fluctuations  on  the  concentration  in  a 
THM  system 
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Pig.  8.  Change  of  temperature  profile  in  the  solvent  zone  for  the  case  of 
a  sudden  temperature  drop 
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Pig.  S.  Heat  flow  and  isotherms  in  THM  for  the  case  of  equal  thermal 

conductivity  in  the  solid  and  the  liquid  [(a)  heat  source  is  shorter 
than  heat  sink,  (b)  heat  source  is  longer  than  heat  sink,  (c)  heat 
source  and  sink  are  equal  in  length] 
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The  volume  of  the  liquid  zone  is  essentially  independent  of  the  temperature 
distribution;  it  is  hounded  by  isothermal  surfaces  at  the  liquidus  temperatures  for 
the  composition  of  the  liquid  (this  would  be  strictly  correct  for  the  static  case  when 
the  zone  is  allowed  to  equilibrate,  but  not  for  the  moving  zone  case) .  With  any  given 
temperature  distribution  (within  a  range  that  should  be  examined) ,  there  will  always 
be  two  isotherms  (equal  or  differing  by  the  AT  required  for  the  growth  rate)  that 
enclose  die  required  volume  of  liquid.  The  problem  is  to  adjust  the  temperature- 
liquid  volume-speed  relationship  so  that  the  isotherm  at  Tg  has  the  optimum  shape; 
in  general,  for  a  given  temperature  distribution,  the  greater  the  amount  of  liquid,  the 
more  concave  the  interlace  becomes  (that  is,  lower  Interface  coucave  to  the  liquid) . 
Also,  for  a  given  amount  of  liquid,  the  lower  the  temperature,  the  more  convex  the 
interface  becomes. 

The  ideal  interface  shape  would  be  planar.  This  seems  to  be  impossible,  even 
ideally,  because  the  system  is  not  sufficiently  symmetrical.  The  most  undesirable 
situation  seems  to  be  a  convex  Interface,  because  the  perimeter  would  be  starved  of 
cadmium,  lower  its  liquidus  temperature,  and  drop  below  the  Tj  isotherm.  Once 
this  starts,  the  instability  would  develop  further,  eventually  reaching  (and  probably 
overshooting)  the  eutectic  point  and  then  forming  the  tellurium  rich  phase.  The  local 
drop  in  interface  temperature  could  also  lead  to  constitutional  supercooling  and  nu- 
cleation  of  new  crystals. 

The  concave  interface  is  less  dangerous,  because  there  is  less  cooling  at  the 
center  than  at  the  periphery.  But  if  the  density  of  the  liquid  increases  with  its  tellu¬ 
rium  concentration,  then  a  pool  of  tellurium  rich  liquid  would  tend  to  form  (by  drain¬ 
age)  in  the  center,  which  could  again  cause  a  progressive  drop  in  liquidus,  possibly 
resulting  in  second  phase  and  the  nucleation  of  new  crystals.  The  implication  is  that 
the  interface  should  be  maintained  slightly  bur.  not  drastically  concave. 

b.  Interface  shape,  partially  mixed  liquid:  Assume  that  the  liquid 
is  stirred  sufficiently  for  all  of  it  (except  the  two  boundary  layers)  to  be  at  the  same 
temperature.  The  heat  flow  situation  is  now  as  shown  in  Fig.  10.  It  would  seem  that 
the  interface  isotherm  could  be  planar  or  concave  but  not  convex  (convex  would 
imply  that  heat  is  extracted  more  efficiently  from  the  center  than  from  the  outer 
ports,  as  heat  is,  in  this  case,  supplied  uniformly  everywhere).  A  slightly  concave 
interface  might  be  satisfactory,  although  the  inherent  disadvantages  mentioned  above 
might  apply.  The  catastrophic  instability  near  the  surface  would,  however,  be  avoided. 
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Fig.  10.  Heat  flow  and  isotherms  in  a  partially  mixed  liquid 
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4.  Conclusions 

The  major  conclusions  from  the  present  analysis  of  the  process  would 
he  that  firsth  the  temperature  gradient  in  the  solvent  (unless  it  became  very  small) 
would  prevent  constitutional  supercooling  at  the  growth  interface.  The  only  other  con¬ 
dition  leading  to  constitutional  supercooling  would  be  severe  temperature  fluctuations 
over  the  short  term.  Thus,  good  short  term  temperature  control  is  essential  and  may 
lead  to  the  necessity  of  controlling  the  voltage  input  into  the  furnace  control  system. 
Tills  would  avoid  temperature  changes  due  to  fluctuations  of  voltage  in  the  mains  which 
can  be  severe  enough  and  fast  enough  to  override  the  control  system.  Secondly,  some 
protection  against  at  least  short  term  power  failures  must  be  provided  in  any  pilot 
production  setup  and  developmental  fornaces  should  at  least  possess  some  means  of 
indicating  short  term  power  failures.  But  beyond  this,  the  inevitable  major  conclusion 
is  that  the  problems  of  inclusions  and  polycrystallinity  may  be  attributed  to  unsatis¬ 
factory  interface  shape;  because  the  temperature  distribution  is  necessarily  unsym- 
metrical,  it  may  not  be  possible  to  achieve  the  optimum  shape  for  both  interfaces; 
the  lower  one  is,  however,  much  more  critical,  although  severe  departure  of  the 
upper  one  from  a  planar  form  would  have  some  effect  on  the  solute  gradient.  There 
is,  however,  no  instability  in  the  shape  of  the  upper  one.  The  analysis  suggests  that, 
if  the  interface  profile  can  be  observed  and  controlled  (by  sufficiently  close  tempera¬ 
ture  control) ,  the  unstirred  liquid  would  be  preferable,  with  a  slightly  concave  inter¬ 
face,  Since  it  seems  not  possible  to  do  this,  due  to  the  fact  that  the  solvent  zone  is 
stirred  by  the  temperature  gradient  Imposed  on  it,  a  minimally  convex  interface  would 
be  desirable. 

B.  Crystal  Growth* 

1.  THM- Experiments  on  the  growth  of  cinnabar 

The  foregoing  discussion,  then,  highlights  the  fundamental  problem 
of  growing  a-HgS  by  the  THM  method.  Due  to  the  relatively  low  temperature  at  which 
cinnabar  transforms  to  metacinnabar,  a  limitation  on  the  growth  temperature  exists; 
the  temperature  at  the  growth  interface  must  not  exceed  the  transformation  tempera¬ 
ture  of  344  °C.  At  these  low  temperatures,  it  becomes  quite  difficult  to  properly  con¬ 
trol  the  thermal  conditions  in  order  to  maximize  the  temperature  gradient  into  the 
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liquid,  which  avoids  constitutional  supercooling  and  stabilizes  the  growth  process. 
Furthermore,  it  is  equally  difficult  to  achieve  a  sharp  temperature  drop  at  the  growth 
interface  in  order  to  cause  the  required  flat  isotherm  there. 

Thus,  it  becomes  extremely  important  to  design  a  furnace  which  allows  these 
thermal  parameters  to  be  optimized.  Such  a  furnace  design,  which  is  now  used  in  all 
growth  experiments,  is  shown  in  Fig.  11.  The  furnace  contains  an  insert  of  p\roiytic 
graphite  with  the  crystallographic  directions  aligned  such  that  radial  heat  flow  is 
facilitated  and  longitudinal  heat  flow  vastly  decreased.  Pvrographite,  as  used  here, 
possesses  a  thermal  conductivity  of  3.9  which  is  better  than  that  of  copper  in  the 
radial  direction  and  a  thermal  conductivity  of  0.035  vtfiich  is  worse  than  that  of 
alumina  ceramic  in  the  longitudinal  direction.  This  combination  of  properties  allows 
the  highest  possible  temperature  gradient  in  the  liquid  in  conjunction  with  the  highest 
possible  drop  at  the  liquid  solid  interface.  The  furnace  profile  is  shown  in  Fig.  12, 
and  a  summary  of  all  crystal  growth  runs  for  the  last  6  months  of  the  contract  is 
found  in  Table  III. 

The  various  heat  flow  difficulties  are  compounded  by  the  fact  that  an  ideal 
solvent  has  not  been  found  yet,  although  we  find  that  which  we  presently  use 
gives  better  results  than  Na^S^  which  is  described  in  the  literature.*  Nevertheless, 
we  believe  that  these  polysulfidcs  probably  possess  quite  high  viscosities,  which 
limit  the  convection  in  the  THM  system,  coupled  with  probably  small  diffusion  co¬ 
efficients  for  Hg,  both  of  these  properties  hampering  effective  transfer  of  mercury 
for  growth. 

With  these  provlslo  in  mind,  we  can  now  examine  Table  I  which  describes 
the  crystal  growth  runs  during  the  reporting  period.  It  can  readily  be  seen  that  the 
use  of  KjS^  as  a  solvent  at  the  highest  possible  interface  temperatures  was  the  most 
promising.  In  several  cases,  bright  rod  transparent  boules  regrew,  which  con¬ 
tained  fat  needlelike  crystals  of  ~  1-2  mm  in  diameter  and  3-4  mm  length  (Fig.  13) . 

No  evaluation  other  than  optical  transmission  through  the  whole  slice  was 
performed  on  there  ingots,  but  such  measurements  must  be  treated  with  caution 
since  the  boule  did  contain  solvent  Inclusions. 

Dr.  A.N.  Mariano  of  Ketmccott  Copper  Corporation's  Ledgemont  Laboratories 
was  kind  enough  to  evaluate  the  growth  habits  of  some  of  these  crystals.  His  con¬ 
clusions  were  that  in  some  of  the  cases  the  habit  of  the  crystals  grown  was  purely 
rhombohedral  with  the  only  form  present  being  { lOll }.  In  other  cases  of  solution 
grown  crystals,  well  developed  forms  of  both  {lOll }  and  (Olll )  were  found. 
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rig.  II.  Design  of  present  THM  ftimaccs 
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Pig.  12.  Typical  temperature  profile  in  a  THM  furnace  as  shown  in  Fig.  11 
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2.  Other  methods  of  crystal  growth  applied  to  cinnabar 


a.  Vapor  growth:  For  the  vapor  growth,  a  modified  Piper  and 
Ifolich  technique  which  had  been  used  at  Tyco  before  for  the  growth  of  CdS  was  used 
(Fig.  14) .  In  it,  the  end  of  an  ampoule  is  slowly  withdrawn  from  a  fUmace  such  that 
a  solid  crystalline  body  grows  from  the  tip  of  the  ampoule,  slowly  filling  the  whole 
space.  The  process  leads  to  growth  of  a  solid  plug  of  HgS,  S  cm  long  in  about  4  days 
under  the  following  conditions:  hot  side:  ~  600  to  700 °C  and  estimated  growth 
temperature:  ~  300-350 0  C. 

Fig.  15  shows  an  ingot  prepared  in  this  way,  sectioned  along  its  length.  Under 
the  microscope  it  can  be  seen  that  this  ingot  is  in  fact  red  colored  and  thus  is  most 
likely  cinnabar.  Although  microscopically  no  evidence  could  be  found  that  the  plug 
had  cooled  through  the  metacinnabar /cinnabar  transformation,  we  still  have  doubts 
whether  all  of  the  material  grew  initially  as  cinnabar.  So  far,  irrespective  of  the 
temperature  gradients  used,  only  very  polycrystalline  ingots  or  collections  of  rods 
were  grown  (Fig.  16) .  Such  rods  or  longer  needles  grown  on  the  walls  were  also 
evaluated  by  Dr.  Mariano  as  to  their  growth  habits.  Longer  needles  were  found  to  be 
well  formed  prisms  with  (lOTO)  habits  twinned  with  (0001 }  as  the  composition 
plane.  Other  habits  found  in  vapor  growth  experiments  were  short  stubby  crystals 
with  well  developed  ( loll }  and  ( Oils }  forms. 

b.  Hydrothermal  growth:  Three  experiments  have  been 

attempted  using  the  hydrothermal  growth  method  of  crystal  growth  discussed  by 

a  10 

Y.  Toudic  and  R.  Aumont  and  also  S.D.  Scott  and  H.L.  Barnes.  The  procedure 
used  was  as  follows:  The  autoclave  (made  by  Autoclave  Engineers,  Inc.,)  was 
thoroughly  cleaned  and  dried  and  the  Teflon  Insert  fitted  with  a  seed  holder  at  the 
bottom  and  2/3  up  from  the  bottom  (see  Fig.  17) .  The  bottom  seed  was  a  slice  of 
single  crystal  of  CdS.  The  top  was  a  Teflon  cup  which  contained  ~  8  g  of  vapor- 
grown  polycrystalline  HgS  material,  A  fresh  solution  of  0.5M  Na^S  was  prepared 
and  ~  If  was  saturated  with  25  g  of  HgS  powder.  The  saturated  solution  was  care¬ 
fully  poured  into  the  Teflon  container  which  fitted  tightly  with  the  autoclave  so  that 
very  little  space  was  left  for  condensate  to  collect.  The  top  flange  was  bolted  down 
firmly  and  the  heating  tape  which  surrounded  the  stainless  steel  vessel  was  turned 
on.  Approximately  24  hr  later  the  temperature  stabilized  ~  230  °C  and  pressure 
of  406  psla.  The  use  of  Teflon  as  a  container  forbids  any  higher  temperature. 
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Fig,  14,  Schematic  of  a-HgS  vapor  growth  furnace 
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Pig.  18.  Longitudinal  section  through  vapor  grown  polycrystalline  a-HgS  boule 


Pig.  18.  Top  of  vapor  grown  ar-HgS  ingot  showing  a  series  of  rods 
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I  Hiring  the  experiment,  the  solution  is  stirred  slowly  (~  30  rpm)  to  aid  the  nutrient 
diffusion  to  the  seed  site.  Results  for  the  first  two  experiments  have  been  negative 
due  to  the  fact  that  the  autoclave  developed  a  leak  after  3  or  4  days  running  and  had 
to  be  shut  down.  In  the  third  attempt,  however,  some  indication  that  a  CdS  seed  may 
work  was  obtained  since  some  IlgS  did  begin  to  grow  on  the  substrate,  despite  a  leak 
which  eventually  terminated  the  experiment. 

3.  Growth  of  materials  other  than  cinnabar 

Some  experiments  on  the  growth  of  BijSj  were  started.  It  has  been 
variously  reported  that  this  compound  has  a  band  gap  between  1.1  and  1.3  eV  and  can 
be  made  with  a  high  resistance  if  excess  sulfur  is  employed.**  Thus,  we  felt  that 
the  high  molecular  weight  coupled  with  a  high  bandgap  and  an  asymmetrical  structure 
might  make  the  material  worth  looking  at  for  electrooptic  purposes.  Thus  BijSj  ingots 
were  grown  by  the  Bridgman  method  and  by  THM  from  Bi  and  subsequently  annealed 
in  a  sulfur  atmosphere.  Also,  growth  from  Te  was  attempted  with  the  aim  of  filling 
the  S  vacancies.  So  far,  the  results  were  disappointing,  since  no  material  of  optical 
transparency  and  high  resistivity  was  obtained. 

Further  experiments  on  the  growth  of  proustite  were  also  conducted  with  much 
the  same  results  as  in  the  first  half  of  the  contract  which  demonstrated  that  proustite 
can  indeed  by  grown  by  THM.  All  these  runs  are  also  summarized  in  Table  111. 


IV.  CONCLUSIONS 


THM  growth  of  cr-HgS  has  proven  to  be  quite  difficult  although  at  present  it 
is  not  definitely  clear  whether  the  method  can  succeed. 

Vapor  growth  of  the  compound  does  not  seem  to  be  practicable  for  large 
crystals  although  it  is  an  excellent  means  to  provide  solid  boules  of  purified  poly* 
crystalline  material  for  regrowth  by  THM. 

Hydrothermal  growth  has  shown  some  initial  promise  and  will  be  pursued 
further.  Also,  some  gel-growth  experiments  are  planned. 

PbCsClj  and  by  inference  PbCsBr^  can  be  grown  if  needed  by  THM  from 
PbClj  (PbBrj)  melts.  Proustite  can  be  grown  from  AsgSg  rich  compositions  in  the 
AggS  -  AsgSg  system.  It  is  transparent  and,  with  some  minimal  farther  work,  could 
be  obtainable  in  good  quality.  If  indeed  the  fact  of  the  disappearance  of  the  10-p 
absorption  peak  in  such  solution  grown  material  can  be  confirmed,  THM  might  pro¬ 
vide  better  quality  material. 

BlgSg  can  be  grown  by  THM  from  B1  solvents.  It  is  obtained  with  good 
crystalline  quality  in  this  case.  Nevertheless,  it  possesses  high  carrier  concentra¬ 
tions  and  therefore  high  1R  absorption.  Other  solvents  such  as  Te  or  annealing  pro¬ 
cedures  have  not  yet  effected  any  significant  change. 
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